A Krylov subspace algorithm for evaluating the phi-functions appearing in exponential integrators

We develop an algorithm for computing the solution of a large system of linear ordinary differential equations (ODEs) with polynomial inhomogeneity. This is equivalent to computing the action of a certain matrix function on the vector representing the initial condition. The matrix function is a linear combination of the matrix exponential and other functions related to the exponential (the so-called phi-functions). Such computations are the major computational burden in the implementation of exponential integrators, which can solve general ODEs. Our approach is to compute the action of the matrix function by constructing a Krylov subspace using Arnoldi or Lanczos iteration and projecting the function on this subspace. This is combined with time-stepping to prevent the Krylov subspace from growing too large. The algorithm is fully adaptive: it varies both the size of the time steps and the dimension of the Krylov subspace to reach the required accuracy. We implement this algorithm in the Matlab function phipm and we give instructions on how to obtain and use this function. Various numerical experiments show that the phipm function is often significantly more efficient than the state-of-the-art.Comment: 20 pages, 3 colour figures, code available from http://www.maths.leeds.ac.uk/~jitse/software.html . v2: Various changes to improve presentation as suggested by the refere

Space-based research in fundamental physics and quantum technologies

Space-based experiments today can uniquely address important questions related to the fundamental laws of Nature. In particular, high-accuracy physics experiments in space can test relativistic gravity and probe the physics beyond the Standard Model; they can perform direct detection of gravitational waves and are naturally suited for precision investigations in cosmology and astroparticle physics. In addition, atomic physics has recently shown substantial progress in the development of optical clocks and atom interferometers. If placed in space, these instruments could turn into powerful high-resolution quantum sensors greatly benefiting fundamental physics. We discuss the current status of space-based research in fundamental physics, its discovery potential, and its importance for modern science. We offer a set of recommendations to be considered by the upcoming National Academy of Sciences' Decadal Survey in Astronomy and Astrophysics. In our opinion, the Decadal Survey should include space-based research in fundamental physics as one of its focus areas. We recommend establishing an Astronomy and Astrophysics Advisory Committee's interagency ``Fundamental Physics Task Force'' to assess the status of both ground- and space-based efforts in the field, to identify the most important objectives, and to suggest the best ways to organize the work of several federal agencies involved. We also recommend establishing a new NASA-led interagency program in fundamental physics that will consolidate new technologies, prepare key instruments for future space missions, and build a strong scientific and engineering community. Our goal is to expand NASA's science objectives in space by including ``laboratory research in fundamental physics'' as an element in agency's ongoing space research efforts.Comment: a white paper, revtex, 27 pages, updated bibliograph

Observational constraints on Horava-Lifshitz cosmology

We use observational data from Type Ia Supernovae (SNIa), Baryon Acoustic Oscillations (BAO), and Cosmic Microwave Background (CMB), along with requirements of Big Bang Nucleosynthesis (BBN), to constrain the cosmological scenarios governed by Horava-Lifshitz gravity. We consider both the detailed and non-detailed balance versions of the gravitational sector, and we include the matter and radiation sectors. We conclude that the detailed-balance scenario cannot be ruled out from the observational point of view, however the corresponding likelihood contours impose tight constraints on the involved parameters. The scenario beyond detailed balance is compatible with observational data, and we present the corresponding stringent constraints and contour-plots of the parameters. Although this analysis indicates that Horava-Lifshitz cosmology can be compatible with observations, it does not enlighten the discussion about its possible conceptual and theoretical problems.Comment: 11 pages, 6 figures, version published in JCA

Classical tests in brane gravity

The vacuum solutions in brane gravity differ from those in 4D by a number of additional terms and reduce to the familiar Schwarzschild metric at small distances. We study the possible roles that such terms may play in the precession of planetary orbits, bending of light, radar retardation and the anomaly in mean motion of test bodies. Using the available data from Solar System experiments, we determine the range of the free parameters associated with the linear term in the metric. The best results come from the anomalies in the mean motion of planets. Such studies should shed some light on the origin of dark energy via the solar system tests.Comment: 10 pages, no figures, to appear in CQ

Freiburg RNA tools: a central online resource for RNA-focused research and teaching.

The Freiburg RNA tools webserver is a well established online resource for RNA-focused research. It provides a unified user interface and comprehensive result visualization for efficient command line tools. The webserver includes RNA-RNA interaction prediction (IntaRNA, CopraRNA, metaMIR), sRNA homology search (GLASSgo), sequence-structure alignments (LocARNA, MARNA, CARNA, ExpaRNA), CRISPR repeat classification (CRISPRmap), sequence design (antaRNA, INFO-RNA, SECISDesign), structure aberration evaluation of point mutations (RaSE), and RNA/protein-family models visualization (CMV), and other methods. Open education resources offer interactive visualizations of RNA structure and RNA-RNA interaction prediction as well as basic and advanced sequence alignment algorithms. The services are freely available at http://rna.informatik.uni-freiburg.de

Single hadron response measurement and calorimeter jet energy scale uncertainty with the ATLAS detector at the LHC

The uncertainty on the calorimeter energy response to jets of particles is derived for the ATLAS experiment at the Large Hadron Collider (LHC). First, the calorimeter response to single isolated charged hadrons is measured and compared to the Monte Carlo simulation using proton-proton collisions at centre-of-mass energies of sqrt(s) = 900 GeV and 7 TeV collected during 2009 and 2010. Then, using the decay of K_s and Lambda particles, the calorimeter response to specific types of particles (positively and negatively charged pions, protons, and anti-protons) is measured and compared to the Monte Carlo predictions. Finally, the jet energy scale uncertainty is determined by propagating the response uncertainty for single charged and neutral particles to jets. The response uncertainty is 2-5% for central isolated hadrons and 1-3% for the final calorimeter jet energy scale.Comment: 24 pages plus author list (36 pages total), 23 figures, 1 table, submitted to European Physical Journal

Measurements of fiducial and differential cross sections for Higgs boson production in the diphoton decay channel at s√=8 TeV with ATLAS

Measurements of fiducial and differential cross sections are presented for Higgs boson production in proton-proton collisions at a centre-of-mass energy of s√=8 TeV. The analysis is performed in the H → γγ decay channel using 20.3 fb−1 of data recorded by the ATLAS experiment at the CERN Large Hadron Collider. The signal is extracted using a fit to the diphoton invariant mass spectrum assuming that the width of the resonance is much smaller than the experimental resolution. The signal yields are corrected for the effects of detector inefficiency and resolution. The pp → H → γγ fiducial cross section is measured to be 43.2 ±9.4(stat.) − 2.9 + 3.2 (syst.) ±1.2(lumi)fb for a Higgs boson of mass 125.4GeV decaying to two isolated photons that have transverse momentum greater than 35% and 25% of the diphoton invariant mass and each with absolute pseudorapidity less than 2.37. Four additional fiducial cross sections and two cross-section limits are presented in phase space regions that test the theoretical modelling of different Higgs boson production mechanisms, or are sensitive to physics beyond the Standard Model. Differential cross sections are also presented, as a function of variables related to the diphoton kinematics and the jet activity produced in the Higgs boson events. The observed spectra are statistically limited but broadly in line with the theoretical expectations

Measurement of χ c1 and χ c2 production with s√ = 7 TeV pp collisions at ATLAS

The prompt and non-prompt production cross-sections for the χ c1 and χ c2 charmonium states are measured in pp collisions at s√ = 7 TeV with the ATLAS detector at the LHC using 4.5 fb−1 of integrated luminosity. The χ c states are reconstructed through the radiative decay χ c → J/ψγ (with J/ψ → μ + μ −) where photons are reconstructed from γ → e + e − conversions. The production rate of the χ c2 state relative to the χ c1 state is measured for prompt and non-prompt χ c as a function of J/ψ transverse momentum. The prompt χ c cross-sections are combined with existing measurements of prompt J/ψ production to derive the fraction of prompt J/ψ produced in feed-down from χ c decays. The fractions of χ c1 and χ c2 produced in b-hadron decays are also measured